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Multifractal properties of wave functions in a disordered system can be derived from self- 
consistent theory of localization by Vollhardt and Wolfle. A diagrammatic interpretation of 
results allows to obtain all scaling relations used in numerical experiments. The arguments 
are given that the one-loop Wegner result for a space dimension d = 2 + e may appear to 
be exact, so the multifractal spectrum is strictly parabolical. The cr-models are shown to be 
deficient at the four-loop level and the possible reasons of that are discussed. The extremely 
slow convergence to the thermodynamic limit is demonstrated. The open question on the 
relation between multifractality and a spatial dispersion of the diffusion coefficient D(uj,q) is 
resolved in the compromise manner due to ambiguity of the D(oj,q) definition. Comparison 
is made with the extensive numerical material. 


1. Introduction 

In previous papers mmm we have initiated a 
systematic analysis of numerical algorithms used in 
the Anderson transition studies [6]. Suggesting va¬ 
lidity of self-consistent theory of localization by Voll¬ 
hardt and Wolfle [7], we have derived the finite-size 
scaling equations for the minimal Lyapunov expo¬ 
nent PP, the mean conductance [2] and level statis¬ 
tics 0. Comparison with numerical results shows 
m m Ei a [g that on the level of raw data they are 
perfectly compatible with the self-consistent theory, 
while the opposite statements of the original papers 
are related with ambiguity of interpretation. It gives 
a serious support to arguments mm that the Voll¬ 
hardt and Wolfle theory predicts the exact critical 
behavior. 

The present paper deals with the next algorithm 
based on the finite-size scaling for inverse participa¬ 
tion ratios [6], which are defined as 

P q = J d d r |T(r)| 2 «, (1) 

where T (r) is a normalized wave function of an elec¬ 
tron in a finite disordered system having a form of 
the d-dimensional cube with a side L. In the metal¬ 
lic state, the wave function 'L(r) extends along the 
whole system and the normalization condition gives 
|T(r)| 2 ~ L~ d and P q ~ £- d (9-i)_ j n the criti¬ 
cal region, the wave functions acquire multifractal 


properties, so 

(P q ) ~ L~ D ^ q ~ 1} ~ L~ d( - q - 1)+A '> (2) 

and the geometrical dimension d is replaced by a set 
of fractal dimensions D q . According to Wegner [IU] , 
the following result takes place for a space dimension 
d = 2 -(- c 

A q = q(q — l)e + 0(e 4 ), (3) 

so the spectrum of anomalous dimensions A q is 
parabolic in the first e-approximation. 

The fractal dimensions D q determine the behavior 
of certain correlators; in particular, 

<|T(r)| 2 |T(r')| 2 )~|r —r'|"\ (4) 

where 

V = d-D 2 . (5) 

Equation (4) is valid in the critical region L < £, 
where £ is the correlation length. In the metallic 
phase, such behavior persists on the scales |r—r'| <£, 
while the constant limit is reached for |r — r'| > £. 
In the dielectric region, dependence (4) is valid for 
|r— r'| <£ and changes by exponential decreasing for 
|r—r'| >£. Since integration of (4) over r and r' gives 
unity, one can estimate the proportionality constant 
in the right hand side and obtain for P 2 

( L~ d (^/a) v (metal) 

(P 2 ) ~ < L~ d (L/a) v (critical region) , (6) 

1 C~ d (U a ) n (dielectric) 
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where a is an atomic scale. Three results in ( 6 ) 
match at £ ~ L, and a comparison with ( 2 ) leads 
to relation (5). 

It is usually accepted [5] that beyond the critical 
point Eq.2 is replaced by the following relation 

(P q )=L- D ^~^F(L/0 , (7) 

which can be used for investigation of the critical 
behavior of £. Below we show that self-consistent 
theory of localization allows to reproduce results ( 2 - 
7) and obtain all functional relations in the explicit 
form. The calculated scaling functions can be com¬ 
pared with the extensive numerical material. Anal¬ 
ogously to HH [3 HE], it appears that the raw nu¬ 
merical data are perfectly compatible with the Voll- 
hardt and Wolfle theory, while the opposite state¬ 
ments of the corresponding authors are related with 
ambiguity of interpretation and existence of small 
parameters of the Ginzburg number type. 

According to certain authors mm, a spatial 
dispersion of the diffusion coefficient D(ui, q) is also 
related with multifractal properties. The diffusion 
constant Dl of a finite system of size L is determined 
for a given function D(uj. q) by the relation 

D L ^ D(D l /L 2 ,L- 1 ) . (8) 

If the power law dependence in w and q is accepted, 
then it is easy to see that a combination 

D(u, q) ~ uj ri ' /d q d ~ 2 ~ r >' ( 9 ) 

provides the correct behavior Dl ~ L 2 ~ d at the crit¬ 
ical point m for an arbitrary value of the expo¬ 
nent rf. The hypothesis put forward by Chalker m 
suggests an equality ?/ = 77 , supported in mm 
by a detailed numerical analysis. In our opinion, 
these arguments are logically deficient: this fact was 
stressed in [14j . but no constructive alternative was 
suggested. 

On the other hand, attempts to introduce a spatial 
dispersion into the scheme of self-consistent theory 
of localization PB CEB reveal the utmost undesir¬ 
ability of this modification. In absence of a spatial 
dispersion, the theory possesses a lot of merits: 

(a) it provides the Wegner relation s = v(d — 2) 
between critical exponents of conductivity (s) and 
the correlation length (V); 

(b) it gives the correct value of the upper critical 
dimension d c 2 = 4, which is a rigorous consequence 
of the Bogoliubov theorem Hz] on renormalizability 
of (j > 4 theory mi; 


(c) it gives the correct dependence D(u, 0) ~ 
w (d-2)/d a f. cr }ti ca J point, which can be obtained 
by different methods HE EOT and was confirmed 
numerically EH; 

(d) it provides a consistent description of finite 
systems considered as zero-dimensional [3]. 

Appearance of a spatial dispersion immediately 
destroys all properties (a-d) El 0: it hardly can 
be considered as incident, since the Vollhardt and 
Wolfle theory is at least a very successful approxi¬ 
mation. In fact, absence of an essential spatial dis¬ 
persion of D(u, q ) was established by the present 
author [9] in the result of a detailed analysis. 

This contradiction can be resolved in the compro¬ 
mise manner, since a definition of D(u>, q) is am¬ 
biguous and allows the ’’gauge transformation” [9]. 
A spatial dispersion is absent in the ’’natural” gauge 
used in [9], but it arises in other gauges allowing 
the equality 77' = i). Unfortunately, it makes unclear 
what gauge corresponds to the observable diffusion 
coefficient; there are indications that in this case the 
equality 77' = 77 is violated (Sec.6) 

2. Two-point correlator 
2.1. Diagrammatic analysis 

Consider the correlator of two local densities of 
states 

K E+UtE {r,r') = {v e+u {t)v e {t')) = 

Y \ip s {r)\ 2 \ip s ,(r')\ 2 6(E + uj- e s )5(E - e s /)\ 

s,s' / 

( 10 ) 

{ip s (r) and e s are exact eigenfunctions and eigenen- 
ergies for an electron in a random potential), which 
is closely related with correlator (4) and can be ex¬ 
pressed in terms of two-particle Green functions 

K E+u ,E{r,r') = ^2 Re [$ flj4 (r, r, r', r') - 

—^(r.rVV)] • (11) 

Here 

$ flyl (ri,r 2 ,r 3 ,r 4 ) = <Gf +w (ri, r 2 )G^(r 3 , r 4 )) 

( 12 ) 

and $ RR is defined analogously. Practically, the 
diagrammatic technique is applied to the quantity 
$^(q) (Fig.l), which is the Fourier transform of 
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Figure 1 : (a) Relation of the function < h£k;'( ( l) with the irreducible vertex (q); ( 6 ) the ladder diagrams; 
(c) a definition of the ” cooperon”. 


( 12 ) with the three-momenta designations taken into 
account 


$ flA (r 1 )r 2 ,r 3 ,r 4 )=L- 3d £ $£ 2 (q> 

k,k',q 


-e ik-(ri-r 3 )+jk'-(r4—r 2 )+iq-( r i- r 2+r 3 — r 4 )/2 

The quantity contains no diffusion poles and 
its contribution is only essential in the zero order 
over the random potential. The quantity is 

determined by the irreducible four-leg vertex U RA 
(Fig.l,a), which reduces to the ’’cooperon” (Fig.l,c) 
in the metallic phase 122 


^kk'(q) = 


2C/ o7 


—iui + £>o(k + k ') 2 


= C/(k + k'). ( 14 ) 


The full 17-vertex differs from (14) by the replace¬ 
ment of the classical value Dq by the exact diffusion 
coefficient D(uj,q) (see Sec. 6 ); here Uo = W 2 a d , W 
is an amplitude of the random potential, 7 = 7 tUqvf 
is an elastic attenuation, determined by the relation 
7 = — ImS^ in terms of the average Green function^ 


(G?) 


1 

E-e k - ’ 


and vf is the density of states at the Fermi level. In 
particular, the one-cooperon contribution to corre¬ 
lator ( 10 ) has a form 

KeLf^ V ') = ^ ReL ~ 3d E p k(q)- 

k.q.qi 


2C/q7 


—iuj + Dqqi 


P_ 


k+q 


i(q) 


,i(2k-qi)-(r-r') 


( 15 ) 


(where P k (q) = G k+q/ 2 G k- q / 2 ) and is easil y cal " 
culated in the ” pole approximation”, when mo¬ 
menta like qi entering the diffusion denominators 
are neglected in slowly varying functions of type 
P_k+ qi (q). In this approximation, one can easily 
calculate contributions to from the ladder dia¬ 
grams shown in Fig.l,b, which have a qualitatively 
different behavior for even ( 2 n) and odd (2n+l) num¬ 
ber of cooperonsd: 


<jj( 2 ”)(r 1 ,r2,r 3 ,r 4 ) = 2 fc„(r 4 -r 2 ) fc*(r 3 -r 4 )- 


■[[/(n-ra))] 2 ", ( 16 ) 

$( 2 n+ i)( ri , r2 , r3 , r4 ) = k n {r 1 -r 4 ) k n (r 2 -r 3 )+ 

+k n (r 1 -r 3 ) fc„(r 2 -r 4 ) [U(r 1 -r 2 ))] 2n+1 . 


1 Below we omit signs of averaging and accept the energy 
variable to be equal E+qj for functions G R and E for functions 

g a . 


2 To retain symmetry of & RA (ri,r2,r3,r4) relative to per¬ 
mutation of r3 and r4, we have added the contributions of 
diagrams with the reversed lower G-line. 
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Here U{ r) is the reverse Fourier transform of (14), 


where 


U{v)=L~ d Y J 


2 U ol 


—ioj + D 0 q 2 


„d -2 ’ 


( 17 ) 


while k n (r) and k n (r) are short-ranged functions de¬ 
fined as 


k n (r)=L~ d J2{G£G£) n+1 e ikr , 

k 

fc„(r) = L~ d J 2 (G«) n+1 (G k ) n e ikr ( 18 ) 

k 

and decreasing as exp(— r/l) on the mean free path l, 
which has the atomic scale near the Anderson tran¬ 
sition. The function $ i? ‘ 4 (ri, r2, r3, Y4) is exponen¬ 
tially small, if all r* are essentially different, and the 
long-range tails arise only in the case of pairwise 
coinciding arguments: the case 17 = r 2 , r3 = r4 
corresponds to correlator (10), while the case ri = 
r 3 , r 2 = r4 (or equivalently iq = r4, r 2 = r 3 ) corre¬ 
sponds to another correlator 

IC E +u,E{r, r') = (Y 4>s{r)4>s(r')ip S '(r)ip s '(r')- 
\ s,s' 


H(r ) = L~ d J2 
q 


e lq ' r 

q 2 + m 2 


to 2 = (-iu)/D(w) = (, 0 £ . ( 23 ) 

The replacement of Dq by D{u>) extends the formula 
obtained for the metallic phase to the whole range of 
parameters, since it corresponds to the replacement 
of the cooperon lines (Fig.l,b,c) by the [/-vertices 
(Fig.l,a); the ladder diagrams are sufficient, since 
all diagrams has the ladder form in terms of the U 
blocks. In correspondence with | 2 j, we neglect the 
g-dependence of the diffusion coefficient, which is 
inessential in the gauge assumed here (Sec.6). In 
a closed finite system, the diffusion constant has the 
localization behavior D(ui) = where £0 d 

if the correlation length of the corresponding quasi- 
zero-dimensional system (2], so the quantity m 2 in 
( 23 ) is finite. Transition to open systems leads to ap¬ 
pearance of the effective damping 70, which is intro¬ 
duced by a change —iu —> —iio + 70 simultaneously 
in — iu> and in D(u>)\ as a result, the finite diffusion 
constant 7 o£od arises in the static limit and the sign 
of the real part can be omitted in (21) and (22). 


2.2. Insufficiency of the pole approximation 


■S(E + w - e s )S(E - e s /)^ . ( 19 ) 

According to ( 16 ), the long-range tails of correlators 
Ke+ui,e{ r, r') and K.e+w,e{*,*') are determined by 
even and odd orders correspondingly. 

Using values of functions fc„(r) and k n { r) at zero 
(the energy dependence of the density of states i'(e) 
is neglected) 

fcn(0) = , k n ( 0) = ^ 2 n+l 0,71 ’ ( 20 ) 

f°° dx r(l/2)r(n+l/2) 

Un ~ 7-00 {x 2 + 1 )"+! “ r(n + 1) 

one has for the essential contributions to (10) 


Contributions ( 21 ) can be easily summed over n, 


T,< !\ U F TO r dt 

K( r,r ) = — Re / . . =, 

K 7-oo Vl + t 2 \/l + t 2 — u 2 


u = 


2 


n(r-r'), 


( 24 ) 


if a 2 is represented as the double integral (see (20)). 
However, this result is practically useless due to in¬ 
sufficiency of the pole approximation. In order to 
clarify a situation, let estimate a value of (21) for 
r = r'. Using the Ward identity [ 7 ] 


AE k (q) = L~ d Y E 4 k' (q)AGk (q), ( 25 ) 

k' 


A'( 2n )(r,r') 
and to ( 19 ) 
/C (2n+1) (r, r') 


—| Vp Re 

7T Z 


2 

■kvfD{uj) 


n(r-r')) 

( 21 ) 


AG k (q) ee G« +q/2 -G^_ q/2 , AS k (q) ee S* +q/2 -S 
and the relation 

AG k (q) = [— u) + e k + q /2 — e k-q/2 + AS k (q)] P k (q), 


27T 2 


jp Re 


7 TVpD{ul) 


n(r—r')) 


2n+i one has for q = 0 


( 22 ) 


L~ d Y (0)P k ' (0) = 1 + iu/ 2 7 ( 26 ) 
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if ImE^ = —7 is assumed to be independent of k.[f] 
For small q and the vertex t/kk' (q), independent of 
momenta, the following relation takes place [22 


and summation over n gives 


£- rf 5> kk '(q)=Pk(q) 

k' 


27 

—iw + D(co, q)q 2 


L d E UoPv (q) = 1 + iujT - D 0 rq 2 , r = I/27. 

k' 

(27) 

There are serious grounds to expect the analogous 
relation in the general case 

L“ d ^£/ k k'( q )-Pk'(q) = 1 +iu}T-D(u},q)Tq 2 . (28) 

k' 

Indeed, the right hand side of (25), being a function 
of k and q, can be practically specified as a function 
of invariants k 2 , q 2 , k • q and allows an expansion 
over the second and the third of them. In the ab¬ 
sence of the k dependence in the left hand side of 
(25), one can average over directions of k and remove 
the odd orders in k • q. As a result, the right hand 
side of (28) contains only even orders in q. A zero 
order term is specified by (26), while the higher or¬ 
ders can be absorbed by a definition of D(ui, q). The 
slow dependence on the modulus of k is removed by 
estimation at k 2 ss e F . 

The equality r = r' in (11) leads to r 2 = in (13), 
so 'Fkk'^) enters in the form summed over k' and 
the radical simplifications are possible due to (28). 
For example, one has for the diagram with two U 
blocks 

^■ d E $ kk'(q) = 

k' 

= L~ 2d E p k(q)C/kk 1 (q)7 , k 1 (q)t7k 1 k'(q)7 3 k'(q) = 

kik' 

= £"‘ , E i ^(q) C/ kk 1 (q)i\ 1 (q)e <t ’ T - D(t ’' ,)T,a = 

ki 

= P k (q)e 2 [^ T_D(l ^ )T<?2 ] . (29) 

Analogously, in the case of n blocks 

L~ d E $ kk' (q) = P k (q)e n [“ T - D( " ; ’ 9)T92 ] (30) 

k' 

3 The k dependence of has no qualitative significance: 
in particular, it is rigorously absent in the Lloyd model, which 
is quite ordinary from the viewpoint of the Anderson transi¬ 
tion. In the general case, neglecting of the k dependence 
corresponds (in the coordinate representation) to the replace¬ 
ment of short-range contributions by the ^-functional ones 
(see Sec.3). 


iAGk(q) +Q{q) rrn 

-iu + D{u,q)q 2 ' { 1 

The exactly such relation follows from the Bethe- 
Salpeter equation (see formula (63) in [9] ), so the 
introduced function D(u, q) can be identified with 
the diffusion coefficient. 

Substitution of (30) in (10,13) gives for the n-th 
order contribution 

A'(n)( 0 , 0 ) = \KeL~ 2d E p k(q)- 

7 T Z z ' 

k,q 

. e n\iuT-D(u,q)Tq 2 ] = ? (32) 

where the latter equality is valid for u —> 0 in the 
vicinity of the critical point, since D(tu,q) turns to 
zero simultaneously for all q [9] . Introducing dimen¬ 
sionless conductance g = v F DL d ~ 2 [T3] and using 
the relation g ~ (L/£) d_2 valid in the metallic phase 
H9. one has from (21) 

A p2n) (r, 0) ~ v 2 f (£/ r ) 2n(d “ 2) . (33) 

The contributions for different n become comparable 
for r ~ ^ and have the order of v‘ F in correspondence 
with (32). It is clear that result (33) is valid for r >£, 
while the r dependence is saturated for r < The 
latter is a consequence of the delicate cancellations 
in the Ward identity (25): summation over k' in the 
infinite limits removes the pole part of Pkk'(q) due 
to its orthogonality to the function AGk'(q) so 
the pole approximation is completely inapplicable. 
For large |r — r'|, summation over k' is effectively 
restricted by the range |k'| < |r — r'| _1 and the or¬ 
thogonality becomes inessential, restoring validity of 
the pole approximation. 

The whole correlator A'(r, 0) is determined by the 
two cooperon contribution for while summation 
of the series of approximately equal terms is neces¬ 
sary for r < £: the expansion parameter u tends to 
unity for |r — r'| — > 0 and the integral (24) diverges 
due to the logarithmic singularity at u = 1. The spe¬ 
cific form of divergency is determined by the char¬ 
acter of saturation of Ad 2r d(r, 0); for example, the 
power law behavior AT( r, 0) ~ r~ 2 P arises, if 

RP 2n \r 1 0) = v 2 F (l-p n r 2 /e) , r<( (34) 
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with coefficients j3 n ~ n 1 . In fact, the dependence 
A'(r, 0) ~ r~ v is expected from the analogy of (4) 
and (10) (Sec.2.4), so 


f r<i 

K(r,0)-v 2 F ~{ . (35) 

where the zero-order contribution is separated. 

In case of the correlator 1C( r, r'), the equality r 2 = 
r4 is valid in (13) for arbitrary r and r', and the pole 
approximation is deficient from the very beginning; 
one always has the result (32), while (22) has no 
range of applicability^. The use of (31) leads to the 
functional form 


which runs in the metallic phase from the minimal 
value £ m in till infinity, and hence £ is not manifested 
as a significant length scale. Then one can conclude 
from (35) that 

77 = 2e (38) 

in correspondence with the Wegner result (see (2- 
5)). We do not expect that a character of the solu¬ 
tion changes at a scale different from £, so the re¬ 
striction eln(L/r) -C 1 is not essential and the two- 
cooperon behavior persists in the metallic phase for 
arbitrary r: 


I\ (r, r') — Up = Up Re 


—n(r-r') 
ng 


(39) 


/C ( r - r ') = ^2 Rei 

q 


27 TUp 


—iuj + -D(w, q)q 2 


iq-(r-r') 


(36) 

corresponding to the first order contribution (see 
(22)), but the diffusion coefficient is defined in the 
different gauge allowing a spatial dispersion (Sec.6). 


In the localized phase one has m = £ _1 and the 
expansion parameter is u = ln(£/r)/g, so the two- 
cooperon behavior holds for £ exp(— g) < r < L. On 
the other hand, for r < £ we expect the same power 
law, as in the critical region. Therefore, the result 
(39) can be extended to the localized phase. 


2.3. Situation for d = 2 + e. 


2.4. Relation of correlators (10) and (4). 


In the spatial dimension d = 2 + e with e <C 1 one 
has for the expansion parameter in (21,22) 

.. — 2II(r) 1 fiPq*.'. 
nupD upD J q 2 


1/r 



1/L 


1 {L/r) e - 1 
9 e 


(37) 


where g = upDL d is a dimensionless conductance 
and the condition 7n<i _1 is accepted, which is valid 
in the metallic state and the critical region. We have 
in mind that a perturbation theory is constructed for 
open systems, where q = 0 is not allowed value for 
the momentum q and the diffusion constant is finite 
in the static limit [2] . Accepting r satisfying the con¬ 
dition e ln(L/r) <C 1 and having in mind that a value 
of g at the Anderson transition is g c ~ 1/e, one can 
see that the expansion parameter u = In (L/r)/g is 
small for the interval Lexp(— 1/e) < r < L both in 
the metallic and critical region. The limiting value 
(29) is not attained and the two-cooperon expres¬ 
sion is valid for correlator (35). This expression is 
not affected by variation of the correlation length £, 


4 By this reason, the argumentation of Sec.2.3 referring to 
the correlator K(r,r') cannot be applied to /C(r,r'). 


In development of the perturbation theory the sys¬ 
tem is assumed to be open, so its conductance g is 
finite and an expansion over 1/g is possible. Inter¬ 
pretation of expressions like (10) in open systems 
suggests broadening of the d-functions to a width 
r > A, where A = 1 /upL d is a mean level spac¬ 
ing. Then each sum over s and s' contains T /A 
terms, and each 5-function gives a factor 1/T. Sug¬ 
gesting that all terms with s = s' (and correspond¬ 
ingly s Y s> ) have the same statistical properties, 
one has for oj = 0 

K e , e ( r,r') « ^ (\YE{r)\ 2 \ipE{r')\ 2 ) + 

+ ^2 (lV’E( r )| 2 |V’E'( r ')| 2 ) ■ (40) 

Assuming for estimate that ipE{r) = ^(r — R) with 
the permanent envelope ^(r) and a random origin 
R, one can replace averaging over disorder by av¬ 
eraging over R and obtain for the second term in 
(40) 

-^ 2 L ~ 2d J d d Rd d R' |T(r — R)| 2 ^(r'—R')| 2 = u% , 

(41) 

while the first term can be estimated as 

4(A/T)(L/|r-r'|) 4 “- d , (42) 
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if 'l'(r) ~ |r| _a and d/2 <2a < d. In the first 
approximation, the zero order contribution u'/ arises 
from the terms with s ^ s' , while the power law 
behavior corresponding to (4) is determined by the 
terms with s = s'. In fact, such decomposition is 
not rigorous because variations of R and R' are not 
independent, so the second term in (40) contains a 
dependence on r — r' (see Eq.44 below). 

A situation is more transparent in the limit of 
closed systems, when F/A —»• 0. Then, for u = 0, 
the vicinity of energy E contains (with probability 
r/A) one level with a certain number so, so only the 
contribution with s = s' = s o remains in sum (10) 
and the second term in (40) vanishes. For transition 
from open to closed systems one should omit the 
zero-order contribution Up, and then correlator (10) 
can be identified with (4) apart from the constant 
factor. Comparison with (39) gives 

A(\^ e (y)\ 2 \^ e (y')\ 2 ) = L~ 2d Re [L e R(r-r ')] 2 , 

( 43 ) 

where A is determined by the normalization condi¬ 
tion (Sec.4). Using the properties of the diffusion 
propagator II(r) (Sec.5), one can easily show that 
result (43) corresponds to the physical expectations 
on correlator (4): the power law behavior |r — r'\~ v 
taking place for |r—r'| < £ changes for |r—r'| > £ by 
saturation in the metallic phase and by exponential 
decreasing in the localized state. The constant limit 
in the metallic phase is determined by the contribu¬ 
tion of the term with q = 0, which is always present 
in closed systems [2]. It disappears in open systems, 
leading to insignificance of the scale £ in the metallic 
phase, which was discussed above. The terms with 
s = s' are the same in correlators (10) and (19), pro¬ 
viding a ground for a hypothesis [23] on the identical 
behavior of these correlators in the critical region^. 

Consider the case of finite frequencies, w A. 
For r ~ A, the first term is absent in the expression 
of type (40) and comparison with (39) gives 

re 1 2 

—n(r-r') . 

7T Q 

(44) 

The power law behavior of the propagator II(r) per¬ 
sists at scales less than L u , where 


L 2d (\ip E (r)\ 2 \ip E +u,(r')\ 2 ) « 1 + Re 


E and E + u> become statistically independent for 
r>L u \ 23]. 

3. Many-point correlators 

Analogously, one can define the n-point correla¬ 
tors 

Jv(ri,r 2 ,...,r n ) = (v Ei {yi) v E 2 (y- 2 ) ... ^„(r n )) 

( 48 ) 

and relate them with many-particle Green functions, 
e.g. for n = 3 


Jv(ri,r 2 ,r 3 ) = --j-h Im [$ a 47 i (n,n, r 2 , r 2 ,r 3 ,r 3 )+ 

+ (ri r 2 ) + (r 2 o r 3 )] , (47) 

and the correlator is determined by the diagrams 
with three G-lines (Fig.2). Selection of diagrams is 
conveniently made in the coordinate representation, 
where the cooperon vertex (14) has a form 


U c ( ri,r 2 ,r 3 ,r 4 ) = U(ri-r 3 ) <5(ri-r 4 ) <5(r 2 -r 3 ) 

( 48) 

and differs from the full four-leg vertex by the re¬ 
placement of short-range functions like k n (r) by the 
^-functions. Analogously, in the analysis of power 
law tails one can use the (5-functions instead of the 
short-range functions G R ( r) and G A (r). According 
to (48), the coordinates r, and r,- corresponding to 
G-lines, coming up from the left to a cooperon ver¬ 
tex, trade places after passing it (Fig.2), and the 
cooperon line gives a factor t/(r, ; — r^). Long-range 
contributions to correlator (46) are determined by 
diagrams, for which the coordinates of all G-lines re¬ 
turn to the same sequence after passing all cooperon 
vertices. The functions t/(r, — Yj ) can enter only in 
even powers, since the coordinates and Yj should 
be transposed even number of times to restore the 
initial sequence. 

Analogously to Sec.2.4, correlator (46) is related 
to the n-point correlator of wave functions. Ac¬ 
cepting the power law dependence on differences 
Yjj = Yi~ Yj , one can write the most general form of 
such dependence 


Lu, = (i/pu) 1/d , (45) 

while for r > L u it changes by exponential decreasing 
(Sec.4); the wave functions corresponding to energies 

5 In the general case, their behavior is surely different, as 
clear from the estimate of type (41). 


(I’4 , (r 1 )| 2 |'F(r 2 )| 2 ... |\I/(r ra )| 2 ) ~ L~ nd 

L\“" / L \ Pn / L \ 7n / L \ Sn 

ri2j \ri 3 J V r 23 ) v ra-l,ra / 
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(p) 


(d) 


Figure 2: In the n = 3 case, the contributions symmetrical over all r,;j are determined by diagrams (a), 
(b), (c) and the diagrams obtained from them by reflections respective the horizontal and vertical axes, (d) 
The diagrams analogous to (a) exist for arbitrary n. 


+permutations of ] , 


(49) following from (1,2), one has 


where permutations provide a symmetry of the ex¬ 
pression over all r%j . Without loss of generality, one 
can accept 

a n > Pn > In > ■ ■ ■ > $n ■ (50) 

Using (49) and the ’’algebra of multifractality” for¬ 
mulated in [23] . one can derive inequalities for A„. 
If all rij ~ L, then correlator (46) allows decompo¬ 
sition [25] FI 

(|(r x )| 2 > (|4'(r 2 )| 2 ) ... (|'I'(r rl )| 2 ) ~ L~ nd , (51) 

so K n = 0. If all = 0, then divergencies in (49) are 
cut-off on the scale a, and comparing with relation 


A„. — a n + fi n + 7„ + ... + S n . (53) 


Using inequality (50) and taking into account that 
there are n(n —1)/2 terms in the sum, one has A n < 
a n n(n — l)/2 which reduces to 


A n < A2 


n(n — 1) 
2 


(54) 


since the exponent a n is n-independent and coin¬ 
cides with r] = A 2 . Indeed, if r-12 <C L and the rest 
rij ~ L , then (49) givesQ 

(l^ , (r 1 )| 2 | v t'( r 2 )| 2 ) (|H / (r 3 )| 2 ) ... (|4/(r„)| 2 ) ~ 


(|4'(r)| 2n ) - L ~ nd+A » (52) 

6 Since it is not quite evident, we accept the following pro¬ 
cedure. Let introduce the scale L u related with a frequency 
and defined in (45); then the functions II(r*j) are exponen¬ 
tially small for rij >Loj (Sec.4) and only the diagrams without 
cooperon lines survive among diagrams like those in Fig. 2; 
these diagrams correspond to decomposition (51). In case 
L u L, the latter decomposition is valid for rij ~ L and 
remains valid approximately, if the scale L u is increased till 
L. 


~L- nd (L/r 12 f n , (55) 

while it is proportional to according to (4). In 
the main e approximation (see (3)), the sign of equal- 

7 The right hand side of (49) may contain less singular 
terms, determined by exponents a n , (3 n , ..., S n whose sum is 
less than A n . If it occurs that d n > a n , then a n = rj and 
a.n < 77 . Hence, the inequality a n < 77 holds in the general 
case, which is sufficient for validity of (54). 


















ity takes place in (54), and correlator (49) is deter¬ 
mined by the most symmetric configuration 

<|T(r 1 )| 2 |'k(r 2 )| 2 ... V ' 

i<j 

(56) 

The contribution 0(e 4 ) in Eq.3, corresponding to 
the orthogonal ensemble, has a structure —an(n — 
1 )(?i 2 — n + 1 with a > 0 (10] [24], and the following 
inequality follows from (54) for n > 1 

n 2 — n + 1 > 3 , (57) 

which is satisfied for n = 2, 3, 4, .... 

The present paper deals with the usual disor¬ 
dered systems like electrons in a random potential, 
which correspond to the Dyson orthogonal ensem¬ 
ble. However, inequality (54) is not related with 
self-consistent theory and has a general character. 
It reduces to an equality for the parabolic spectrum 
A q = /3q(q — 1) with arbitrary /3, leading to the 
symmetric form (56) with q = 2/3 for the n-point 
correlator. In particular, it is actual in the first 
e-approximation for the unitary ensemble, where 

mm 

A n = n(n-l)(e/2) 1/2 +(3/8)n 2 (n-l) 2 C(3)e 2 , (58) 

and in the regime of the quantum Hall effect (see 
below). Substitution of (58) into (54) gives the in¬ 
equality for n > 1 


n(n — 1) < 2 , (59) 

which is violated for n = 3, 4, ... Hence, the results 
obtained in the tr-models reveal deficiency on the 
four-loop level. A possible reason of that is discussed 
in Sec.7. 

Expression (56) allows a diagrammatic interpreta¬ 
tion. For small e one can neglect non-symmetrical 
terms, though a mechanism of their compensation 
is not quite clear. In case n = 3, the lowest order 
symmetrical contribution arises from the diagrams 
in Fig.2,a-c: 

/ J_d—1 

Jr(ri,r 2 ,r 3 ) = const y% - 

V 9 

■n(r 12 ) 2 n(r 13 ) 2 n(r 23 ) 2 . (60) 

The diagrams analogous to that of Fig.2,a exist for 
arbitrary n, as illustrated in Fig.2,d for n = 4: the 
first (n — 1) cooperons provide a cyclic permutation 



of ri, r 2 ,..., r„, which should be repeated n times, 
in order to restore the initial configuration. Adding 
the zero-order term, one has 

K( ri,r 2 ,...,r„) - ^ = (61) 

/ Td -2 \ 

= const Vp ( -J JJ n (nj ) 2 . 

\ 9 J i<:j 

In the limit of closed systems the term Vp disappears 
(Sec.2.4), and Eq.61 gives the main symmetrical con¬ 
tribution in the metallic region for r t j >£, which can 
be extended to arbitrary ry,- analogously to Sec.2.3. 

Formally, expression (61) is obtained for d = 2 +e 
with small e, but in fact its validity is related with 
two qualitative moments: 

(i) insignificance of the correlation length £ in the 
metallic phase as a characteristic length scale; 

(ii) realization of the maximally symmetric form 
(56) for the n-point correlator. 

These properties can be approximate and valid only 
for small e. However, their qualitative character al¬ 
lows to assume that they persist in the general case, 
as supported by a diagrammatic interpretation of 
results. In such a case, the multifractal spectrum is 
determined by the relation A n = n(n — l)e and ap¬ 
pears to be strictly parabolical. Below we use Eq.61 
in the 3 D case, considering it as extrapolation from 
small £ to e ~ 1, but having in mind that it can be 
exact. 

The simplest arguments do not allow to reject this 
hypothesis. A reference to the e-expansion is un¬ 
founded, since cr-models are deficient on the four- 
loop level (Sec.7). Numerical data are not reliable 
due to extremely slow convergence to the thermody¬ 
namic limit (Sec.5). On the other hand, the follow¬ 
ing arguments can be given in favor of the hypothe¬ 
sis. 

(a) The result q = 2e looks plausible, since the 
condition q > d is fulfilled for d > 4; then it follows 
from (56) that the normalization integral is deter¬ 
mined by the atomic scale for all actual correlators. 
It agrees with estimates by the optimal fluctuation 
method and instanton calculations (25j . which pre¬ 
dict localization of wave functions for d > 4 at the 
atomic scale even in the critical region. 

(b) A surprising accuracy of the Wegner one-loop 
result (3) in application to the d = 3 and d = 4 cases 
was reported in a lot of numerical experiments [261 
[27] EH [29]. For example, a position of the maximum 
for the singular spectrum f(a ) (which is ao = d+e in 
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the one-loop approximation) was estimated as «o = 
4.03 ± 0.05 !l|, a 0 = 4.048 ± 0.003 [27] for d = 3 
and ao = 6.5 ± 0.2 |26; for d = 4. A parabolic form 
of the spectrum is confirmed on the level of 10 % m 
1251 125] , which should be considered as satisfactory 
(Sec.5). 

(c) In the regime of the integer quantum Hall ef¬ 
fect, the spectrum is parabolic on the level of 10~ 3 
m, and there are theoretical arguments in favor of 
strict parabolicity [31. [52] [33] (confirming the prop¬ 
erty (ii)) based on the relation with the conformal 
field theory. 

(d) Validity of the Vollhardt and Wolfle the¬ 
ory is directly related with the property (i). In¬ 
deed, it is known from finite size scaling that g = 
g c + const (L/£) 1 / 1 ' in the critical region [ 2 ] [5], while 
g ~ (L/^) d ~ 2 in the metallic phase [19j[l; it gives the 
relation = d— 2 [7], if £ is not a significant length 
scale. Thereby, in the framework of self-consistent 
theory the property (i) is naturally considered as ex¬ 
act. 

(e) Application of the ” algebra of multifractality” 
to correlators of the more general form than (49) 
leads to the statement on strict parabolicity of the 
multifractal spectrum [311 • Therefore, the symmet¬ 
ric form (56) is exact, while deficiency of er-models 
takes place not only for unitary, but also for the or¬ 
thogonal ensemble. 


side over rq, ..., r„ gives unity: 

A = L -dn+n(n-l)e j ^ J n ( p ._ rj )2 _ 


i<j 


(63) 

Integration is easily performed in case n = 2, giving 
A as a regular function of z = L /£ozq 


A = A(z) = £ 


l 2 


Z 1 + (27rs) 2 


1 / 2 4 , 2 < 1 

d ~ 4 z » 1 


CdZ 


(64) 


where Cd = nKd(l — d/ 2 )/ 2 sin( 7 rd/ 2 ), Kd = 
[2 d-1 7r d/,2 r(d/2)] 1 , and s = (si,..., Sd) is a vector 
with integer components Si = 0, ± 1 , ±2 .... Accord¬ 
ing to [ 2 ], the quantity z is a function of the ratio 
L/£ determined by the equation 


±c d (L/0 d ~ 2 = H(z) , 


(65) 


where Cd = nKd/\2 sin( 7 rd/ 2 )| and H{z) is a function 
introduced in [2 with the asymptotics l/z 2 for 2 -C 1 
and —CdZ d ~ 2 for z> 1. Setting r = r' in (43) and 
substituting to (1), one has for P 2 

(P 2 ) = A " 1 L~ d ( L/af e , ( 66 ) 


4. Scaling for inverse participation ratios in accordance with (7). Using (64),(65), we have 


We have established that the critical behavior of 
correlator (56) is reproduced by the diagrammatic 
contribution (61). The latter has a more wide range 
of applicability and allows to extend the results be¬ 
yond the critical region. In the limit of closed sys¬ 
tems, one has from (61) analogously to Sec.2.4 

(|T(r 1 )| 2 |4>(r 2 )| 2 ...|vl/(r rl )| 2 )= (62) 

= A~ 1 L~ dn l/d ”- 1 ) 6 J|n(r w ) 2 , 

i<j 

where the parameter A is determined by the normal¬ 
ization condition, since integration of the left hand 

8 According to one-parameter scaling theory m, 9 = 
F(L/£) where the function F(x) has a behavior x d 2 in the 
metallic phase, in order to provide the relation g oc L d ~ 2 . 
Due to the dependence £ ~ t~ u (t is a distance to the transi¬ 
tion), one can consider g as a function of the argument tL 1 / u , 
which allows the regular expansion in r due to the absence 
of phase transitions in finite systems; the first order in r is 
sufficient in the critical region. 


( L~ d (£/a) 2e (metal) 

(P 2 ) ~ < L~ d (L/a) 2e (critical region) (67) 
l £ -d (£/ a ) 2e (dielectric) 

in agreement with ( 6 ). 

In case of arbitrary n, one can obtain from (63) 
that A(z) ~ z ~ 2n ^ n ~ 1 ^ in the metallic region and 
A(z) ~ z -d(n-l)+n(n-l)e j n l oca li zec 4 one . The 

first result is a consequence of the fact that propaga¬ 
tor (23) in the region m <C A -1 is dominated by the 
term with q = 0 and practically constant. To obtain 
the second result, one changes from variables rq to 
variables iq and r' = iq — iq (* > 2 ) and exploits the 
iq-independence of the integrand and its localization 
in the |r'| < ^ 0 z> region. Then using (65) one has 

f ~ (L/^) n< ' n ~ 1 ' )£ (metal) 

A = < A c ± B(L/^) d ~ 2 (critical region) . 
[ ~ {L/0~ d{n ~ i)+" ( ”-i )e (dielectric) 

(68) 
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Setting Tij = 0 in (62) one gets in analogy with (66) 


3 


(P„) ~ A~ 1 L~ d I"" 1 ) (L/af (n_1)e , (69) 


and 

ln(P„) = —D n (n - 1) In (L/a) + const + F n (L/£) , 

, (70) 

where the constant is chosen from the condition 
F„(0) = 0 and 


F n (x) = - \n{A/A c ) = 


—n(n — l)elnx 
±B n x d ~ 2 
D n {n — 1) In a: 


(metal) 

(critical region) 
(dielectric) 


(71) 


Evaluation of F n (x) for arbitrary x can be made 
rewriting (63) in a form of the multiple sum over 
momenta. Unfortunately, such expression needs te¬ 
dious numerical calculations for large n and does 
not provide analytic continuation to non-integer n. 
To avoid these problems, we note that the result 
A ~ z —2n(n—\) f or ^he me tallic phase remains valid 
in the critical region by the order of magnitude; 
therefore, in these regions F n (x) differs from F 2 (x) 
by a factor n(n — l)/2, while in the deep of the lo¬ 
calized phase F n (x) = (n — l)(D n /D 2 )F 2 (x). The 
simplest interpolation form ensuring such properties 
is as follows 


J C+F 2 (ax) (upper branch) . 

( C-F 2 (x) (lower branch) ’ 


i.e. two branches of F n (x) have the same behavior 
as two branches of F 2 {x) and differ only by a scale 
transformation. The coefficients C+ and C_ provide 
the correct asymptotic behavior (71) for large x, and 
the parameter a ensures symmetry of two branches 
for small x: 


^ D n (n — 1) ^ 

C+ " Jh ’ C - 


n(n — 1) 
2 



(73) 

If the function F 2 (x) is calculated (Fig.3), one can 
compare (72) with results by Brndiar and Markos 
for n = 5 [35] in three dimensions (Fig.4,a). Due to 
the presence of large parameter n(n — 1) = 20, all 
numerical data lie in the critical region 1, where 
the dependence F n (x ) is practically linear in accor¬ 
dance with v = 1 in the Vollhardt and Wolfle theory. 
Linearity of dependencies in Fig.4,a is also evident, 


2 

1 

0 

-1 

-2 



Figure 3: Scaling function F n (L/£) for n = 2. 


and their matching with the theoretical scaling curve 
offers no problem (Fig.4,b).@ 

The opposite situation takes place for numerical 
data by Rodriguez et al [27; for n = 1/2 (Fig.5,a). 
In this paper the relation was accepted 

<P„) ~ L -d(n-l)+A n ( 74 ) 


for the whole range of parameters, implying that 
A n = A„ at the critical point; then a comparison 
with (69) gives 


A 


n 


A n + 


Fn(L/Q 
In (L/a) 


(75) 


In addition, roughening was made at the length scale 
l, which should be used instead of a in (75). If 
A = l/L is fixed, then (75) contains only a depen¬ 
dence on L/£, which is determined by the function 
F n (x). Numerical data for n = 1/2, A = 0.1 [27 
are shown in Fig.5,a: the stationary limit is reached 
for a value W c = 16.6, which gives the estimate of 
the critical point. Accepting A n (W, L) — A n (W c , L) 
as deviation from the critical behavior, one can 
match all numerical data with the theoretical scal¬ 
ing curve by a change of scale along the horizontal 
axis (Fig.5,b). Due to existence of the small param¬ 
eter n(n— l)/ln(l/A) = 0.11, the main body of data 
corresponds to large values of x = L/£, so the lower 
branchial is determined by its logarithmic asymp¬ 
totics, while the upper branch remains in the linear 

9 Details of the scaling procedure were discussed in Sec.3 
of 0. 

10 Due to the negative value of the factor (n — 1) the upper 
and lower branch trade places, and to restore their natural 
disposition we consider the quantity —A n . A definition of 
A n in |27| was accepted with the opposite sign, and Fig.5,a 
directly corresponds to Fig.6,c in na¬ 


il 









Figure 4: Numerical data by Brndiar and Markos for n = 5 extracted from Fig.2 in [55] . and their 
comparison with the theoretical scaling dependence. The empirical values P >2 = 1.28 and P 5 = 0.96 [33 
were used. A difference between (lnP„) and ln(P„) was neglected. 


0.5 

0.4 

0.3 

0.2 

0.1 

20 40 60 80 100 L 0 2 4 6 8 10 

(a) ( b) 

Figure 5: Numerical data by Rodriguez et al for n = 1/2 extracted from Fig. 6 ,c of [27], and their 
comparison with the theoretical scaling dependence (72). Fig.5 ,b contains points corresponding to L = 
20, 40, 60.... 
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regime due to a small value of a. It explains why 
dependencies for W > W c are practically linear (see 
inset in Fig.5,b), while a tendency to saturation is 
evident for W < W c . Small deviations in Fig.5,b are 
probably related with inaccuracy of the interpolation 
form (72). The evident linearity of dependencies at 
small L corresponds to a value v = 1 of the Volllrardt 
and Wolfle theory, while a statement of m on the 
result v = 1.590 obtained with ’’unprecedented pre¬ 
cision” looks rather strange^- 

For frequencies w A, the following equations 
are valid for g and z = L/£ qd [3] 


value, and one can use the periodical boundary con¬ 
ditions, accepting q = 2ns/L, where s is a vector 
with integer components Sj. For r ^ 0, the sum over 
q is convergent and no cut-off is necessary at large 
momenta. Then one can obtainF^I 

n(r) =]rn 0 (r + sT), (78) 

s 

where Ilo(r) is a continual version of (23) 

d d q e iq ' r 
(27r) d q 2 + m 2 


9L = 


P_ 

*2 


±Cd(L/£) d 2 = - CdZ d 2 , 


(76) = 


_ - _(_LV 

(47r)d/2 \2 m) 




// = 1 - d/2 (79) 


where p = (—iw)/A. At the critical point one has 
^ = oo and z ~ p 1 / d , so £od coincides with the 
scale L u introduced in (45). Therefore, m~ l ~ L w 
and the propagator II(r) falls exponentially on the 
scale L w <C L providing statistical independence of 
|t/tE(r)| 2 and |f/>E+ w (r)| 2 for r and fulfilment of 
the normalization condition in (44) apart to small 
deviations. At the critical point, Eqs.76 give g ~ 
(w/A)^ _2 ^ d in correspondence with 18, [19] [20] ; 
substituting to (44) and setting r = r', one has for 
small frequencies 

{\^E{jc)\ 2 \ip E+UJ {T)\ 2 ) ~ L~ 2d oc u}- v/d . 

(77) 

Numerical verification of such scaling was carried 
out in papers m [38] and was interpreted as con¬ 
firmation of Chalker’s hypothesis m on a spatial 
dispersion of D(u, q). We see that this result can be 
obtained without any reference to the q dependence. 

5. Convergence to the thermodynamic limit 

According to (61), all actual correlators are deter¬ 
mined by the diffusion propagator II(r) defined in 
(23), which should be estimated for closed systems 
(Sec.2.4). The latter contain q = 0 as an allowed 

11 The paper m exploits the treatment procedure devel¬ 
oped in m, which was already criticized m- It involves 
many-parameter nonlinear fitting, which leads to the huge 
number of the \ 2 minima and allows to obtain any desired 
value of v in a rather wide interval. A ” desired” value 
v = 1.590 was chosen from the correspondence with preceding 
papers (occurrence of such values was discussed in DO). while 
its ’’unprecedented precision” corresponds to fluctuations in 
the single x 2 minimum and has no relation to actuality. Anal¬ 
ogous arguments are valid in respect to accuracy of ao (Sec.3) 
and D 2 (Sec.5). 


(Kn(x) is the Mac-Donald function) with the asymp¬ 
totic behavior for d > 2: 


n 0 (r) = < 


r(d/ 2 -l) d _ 2 

(4n) d / 2 

(tt/2) 1/2 m d ~ 


2 

mr 


{2n) d / 2 (mr)( d_1 i/ 2 


- , mr <C 1 


mr 1 


( 8 °) 

To provide a finite value for r = 0 we accept the 
spherical cut-off |q| < A, so 


n 0 (o) 


K d A d ~ 2 
d -2 


(81) 


and the growth at r —> 0 in (80) saturates for r<A -1 . 
According to (78), II(r) is a sum of spherically sym¬ 
metric functions originated in the centers of cubi¬ 
cal blocks with side L. This fact, along with the 
cut-off |q| < A, leads to distortion of dependencies 
(80) specific for infinite systems: it is manifested in 
anisotropy over directions of r and in oscillations in¬ 
duced by the cut-off. As a result, the exponent in 
dependence II(r) oc |r| _ “ is determined in a finite 
system with unavoidable restricted accuracy. 

According to [2] , the relation mL = zq takes place 
at the critical point, where zq is a root of the func¬ 
tion H(z) (see(65)). The parameter zq is not univer¬ 
sal, but depends on the details of cut-off, and hence 
on the specific model (zo ~ 2 for a spherical cut¬ 
off). Fig.6,a illustrates the results for II(r) in the 
3 D system of size L = 20 for A = n and different zq. 
One can see that satisfactory power law dependen¬ 
cies II(r) ~ |r| -a with a = 0.27 — 1.54 take place for 
Zo = 1 — 4, and their quality does not allow to dis- 

12 It follows from the a-representation and the Poisson sum¬ 
mation formula (see Appendix in 0). 
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tinguish a theoretical value a = 1. If L is increased 
to 100 (Fig.6,b) the range of a becomes somewhat 
narrower for the same conditions, but remains rather 
wide (0.36—1.46). Surprisingly, the picture does not 
change essentially even in the thermodynamic limit 
L —> oo (Fig.6,c), if a value of |r| is a finite fraction 
of L. Indeed, setting r = y L, mL = zo , one has from 
(78) (the integration limits are shown for a modulus 
of q) 


A L 


H(r )=L 2 ~ d J2 


d d q 


p *q-(y+s) 


(2Tr) d q 2 + z { 


2 > 


(82) 


so that A L —> oo for large L , and the cut-off is re¬ 
moved effectively; then L enters only in the common 
factor and does not affect the y dependencies. The 
theoretical value a = 1 should appear in the limit 
\r\/L —> 0 independently of z 0 i but in fact this limit 
is unattainable even for the maximal sizes L ~ 100 
reached in the present time [53] [55] , One can see 
from Fig.6,c that for |r| comparable with L the ex¬ 
ponent a is determined by the Zq value specific for 
a given model and has the scattering 0.55 — 1.43 for 
zo = 1 — 4. The satisfactory power law dependen¬ 
cies are observed for |r| > 0.05T, i.e. for |r| > 5 at 
L = 100, while the scales |r| < 5 always fall out a 
scaling picture due to influence of the cut-off. Ac¬ 
cording to (61), II(r) determines a behavior of the 
n-point correlators, and it does not look improbable 
if numerical estimations of fractal dimensions may 
by tens of percents deviate from the true ones. In 
particular, it is hardly possible to make the state¬ 
ments of principle, relying on deviations in the third 
digit [30]. 


Figure 6: Behavior of the diffusion propagator II(r) 
with different Zo for L = 20 (a) L = 100 (b) and 
L = oo (c). The origin at the vertical axis is chosen 
arbitrary. 


T a b l e 1 

Estimations of 77 and D -2 = d — 77 for the 3 D 
systems of different sizqlj L. 


13 In respect of the last result see Footnote 4 in m 
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L 

1 

d 2 

Paper 

10 

1.4 ±0.1 

1.6 ±0.1 

m 

16 

1.3 ±0.3 

1.7 ±0.3 

m 

20 

1.67 ±0.02 

1.33 ±0.02 

m 

40 

1.3 ±0.2 

1.7 ±0.2 

m 

40 

1.5 ±0.3 

1.5 ±0.3 

m 

47 

1.32 

1.68 

m 

48 

1.48 ±0.11 

1.52 ±0.11 

m 

60 

1.38 ±0.18 

1.62 ±0.18 

m 

80 

1.70 ±0.05 

1.30 ±0.05 

m 

120 

1.76 ±0.03 

1.24 ±0.03 

m 

240 

1.76 ±0.07 

1.24 ±0.07 

m 


According to the relation r] = 2a (see (58), one 
can expect the scattering r) = 1.1 — 2.8. Table 1 
gives estimations of r) and D 2 = d — 7] obtained for 
d = 3 by different authors. One can see their large 
scattering and a systematic drift with increasing of 
L. The last estimate rj = 1.76 ± 0.07 has only the 
10% deviation from the value rj = 2 following from 
the first e-approximation, which has a chance to be 
exact according to above arguments. The observed 
deviations from the parabolical spectrum are also on 
the level of 10% [27] M [29]. 

6. On a spatial dispersion of the diffusion 
coefficient 


It is clear from the above discussion that all the 
picture related with multifractality of wave func¬ 
tions can be obtained without any reference to the 
g-dependence of the diffusion coefficient D(uj,q). 
At first glance, it indicates a complete failure of 
Chalker’s hypothesis El- In fact, a situation is more 
complicated due to ambiguity of the D(w,g) defini¬ 
tion. 

The arguments of Hung are based on the relation 


/C(g) = 


i'F D{u, q)q 2 


77 oj 2 ± [D(u],q)q 2 Y 


(83) 


for the Fourier transform of correlator (19) and an 
assumption on the similar behavior of correlators 
JC(r,r') and AT(r,r') in the critical region [23]; then 
starting from A'(r, 0) ~ r~ v one has /C(g) ~ K{q) ~ 
q~ d +v an d D(w,q) ~ q d ~ 2 ~ v . In fact, the correct 
relation has a form (see below) 


and reduces to (83) under assumption of the real 
diffusion constant and irrelevance of the regular 
part <j>reg(q)■ The identical behavior of /C(r,r') and 
A'(r, r') can be ensured in the limit of closed sys¬ 
tems (Sec.2.4) when D(oj, q) oc (— iw) and the pole 
term in (84) gives no contribution in the main order 
of oj. In the general case, complex-valuedness of the 
diffusion coefficient does not allow to draw reliable 
conclusions relative D(w,g) from the given behavior 
of /C(g). 

According to [5], the use of the Kubo formulas 
allows to relate the Fourier transform of (12) for ri = 
r 3 , r 2 = r 4 with the observable diffusion coefficient 


0(g) 


2'ki'f 

—ico ± D(u >, q)q 2 


± 0reg (g) • 


(85) 


Substitution of (85) into the expression for /C(r,r') 
analogous to (11) gives Eq.84 where the regular part 
is somewhat different from (85) due to a contribution 
of Decomposition into the pole and regular- 

parts is not unique and allows the ’’gauge transfor¬ 
mation” [S] 

0re«?(g) = firegid) ~ 2? TV F C{q), 


D(uj,q)q 2 


D(w, q)q 2 + iuC(q)[—iu + D(u, q)q 2 ] 
1 + C(q)[—ioj + D(ui, q)q 2 ] 


( 86 ) 

where C(q) = 0(q 2 ) for small q. Another represen¬ 
tation for (j>(q) follows from the spectral properties 
of the quantum collision operator [9j: if A s (g) are its 
eigenvalues, then 


0(g) 


^o(g) 2 Als(g ) 2 

—w + Ao(g) s ^q— W + A s (g)’ 


(87) 


where A(j(g) = -27 nv F + 0(g 2 ), A s (q) 2 = 0{q 2 ). 
The eigenvalue with s = 0 has a behavior Ao(g) ~ g 2 
for small g and one can accept by definition 

Ao(g) = -iD(u,q)q 2 . (88) 


Then (87) coincides with (85) where the regular part 
behaves as g 2 for small g and can be excluded by the 
gauge transformation. The gauge (88) will be re¬ 
ferred as ” natural”; it was exploited in [9] and found 
to be free of an essential spatial dispersion. An¬ 
other distinguished gauge is defined by the condition 
0reg(g) = 0; it is actual in the localized phase, where 
D(cu,q) = (— ioj)d(q ) and the following relation 


£(«) = ^2 Re 


2ttu f 


-ioj + D(uj , q)q 2 


rreg 


(g) 


(84) 


1 

1 + d(q)q 2 


A(q) 


J dre- zq r A(r) 
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“4( r ) = ^{J2 \Mr)\ 2 \^m 2 S(E - e s )J (89) 

can be obtained from the Berezinskii-Gor’kov crite¬ 
rion jjjj. The term with s = s' in correlator (19) 
gives contribution S(u>) in the localized phase, trans¬ 
forming to the singularity l/to in the quantity 4>(q), 
which can be identified with the diffusion pole in 
(85). It is of principle importance to gather all con¬ 
tributions ~ 1 /lo which may be contained in (f> reg ( q ) 
and to include them in the pole term. They are cer¬ 
tainly present in ip reg (q) for gauge ( 88 ), since there 
are terms with A s (q) ~ w in the sum of (87) [9]. 
Such contributions are surely included, if the gauge 
with 4> reg (q) = 0 is chosen, and just this gauge is 
implied in (89). Comparison of (89) with (4) shows 
that A(r) ~ r~ v for r < £ and d(q) ~ g d_2_7 ) for 
q^ZZ. 1 in correspondence with the Clialker hypoth¬ 
esis! 14 l. while d(q) = const = £ 2 for q < £ _1 and 
A{r) ~ exp(— r/£) for r > £. If the gauge with 
d{q) = const is used, then the contribution q~ d+v /co 
is contained in (f> reg (q ). 

The ’’natural” gauge ( 88 ) was used in the analysis 
of [9j and just that very definition of D(co,q) is im¬ 
plied in the vertex U RA . If the pole approximation is 
used in Eq.12 of (9], then one can set k' = k in the 
function F( k, k', q) and obtain F(k, — k, q) = 2 f 7 o 7 
from Eq.65 of this paper, if ImE^ = —7 is accepted 
to be k-independent and Uq is defined by the relation 
7 = ttUqI'f• Then the pole term of the vertex U RA 
can be obtained from the cooperon contribution (14) 
by substitution of D{w, q) for D 0 and neglecting the 
q dependence. The precisely such form of the vertex 
was used in the above considerations, and its validity 
is confirmed by successful reproducing of multifrac¬ 
tal properties. It gives an essential support to the 
conclusions of 0 . 

The paper jTTJ contains inaccuracy consisting in a 
wrong interpretation of the ’’regularity” of <f> r eg{q)- 
It was suggested that the Anderson transition is 
completely determined by diffusion poles, while the 
function ip reg (q) does not feel the transition and con¬ 
tains no information on the correlation length £. 
Then the quantity C(q) in ( 86 ) relating two such 
regular functions is also ^-independent, and absence 
of the anomalous dispersion (related with the scale 
£) in one gauge leads to its absence in another gauge. 

14 The results for l)(uj. q) in the localized phase and the 
critical region are matched at § ~ ; thus at the critical 

point one has D(oj, q) ~ uj /ijt ~ uk d 2 )/ d for < 1 and 
D(uj, q) ~ uLl(qL u ,) d - 2 - rl ~ ( WV _2_ ’ 7 for qL u > 1. 


In fact, as we see, information on £ is unavoidably 
present either in D(ui,q), or in <p reg (q). 

The latter leads to the disappointing conclusion, 
that nothing can be said on the gauge corresponding 
to the observable diffusion constant, so the exponent 
rf in (9) is in general different from 77 . The most de¬ 
tailed verification of the relation rf = 77 was carried 
out in the paper [T2] . In fact, a value rj' = 1.20±0.15 
was found for the exponent t/ and two estimates 
(77 = 1.3 ± 0.2 and 77 = 1.5 ± 0.3) were obtained for 
the exponent 77 : according to the authors, it was suf¬ 
ficient to establish the equality 77 ' = 77 . Lately, the 
estimate for 77 has drifted to 1.76 (see Table 1), while 
no fresh data for 77 ' is known to us. One can see from 
Fig.7 that the raw data of |T3] for the autocorrelation 
dependence t~^ d ~ v ^ d are perfectly compatible with 
a value 77 ' = 1 corresponding to absence of the spa¬ 
tial dispersion. The physical experiment on spread¬ 
ing of the wave packet [JSj is also well-described by 
self-consistent theory of localization. 

7. On deficiency of er-models 

In Sec.3 we have established deficiency of the a- 
model approach beyond the first e-approximation. 
This situation is not unexpected: derivation of er- 
models is justified only for small e, and the question 
on their exact correspondence with the initial disor¬ 
dered systems always remained open. In particular, 
strong doubts arouse in relation with the upper crit¬ 
ical dimension [5]. Here we present arguments, why 
deficiency of er-models for the orthogonal ensemble 
arises just on the four-loop level. 

Following arguments of Sec.3, we can assume, that 
the two-cooperon form of correlator K (r, r') is exact. 
Then absence of the spatial dispersion of D(w, q) cor¬ 
responds to the exact equality 77 = 2 e, while viola¬ 
tion of this equality corresponds to appearance of 
the spatial dispersion. 

Wegner’s result (3) is obtained in the ’’minimal” 
cr-model, where the lowest (second) powers of gra¬ 
dients are retained: it corresponds to neglecting a 
spatial dispersion of D(u,q). In the first three or¬ 
ders in e, the equality 77 = 2 e takes place and the 
approximation is self-consistent. This equality is vi¬ 
olated on the four-loop level and signals on the lack 
of self-consistency. The terms with higher gradients 
should be taken into account, which leads to insta¬ 
bility of the renormalization group due to the ’’gra¬ 
dient catastrophe” [46] . To remove instability, addi¬ 
tional counter-terms should be added; it leads to es¬ 
sential modification of the cr-model Lagrangian and 
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t/*H 


Figure 7: The raw data by Brandes et al ng on spreading of the wave packet; the autocorrelation 
function C(t) ~ t~^ d ~ v ^ d describes the change of the amplitude in the packet center as a function of 
time. Dependence f -2 / 3 corresponds to the absence of a spatial dispersion of the diffusion coefficient, while 
dependence t ~°' 6 is the result indicated in [12j. 


makes unclear a fate of the four-loop contribution. 
It should be stressed that according to the analysis 
of the paper [5] the spatial dispersion is determined 
by the atomic scale. It is inessential in the practical 
sense but its existence is a matter of principle due to 
the infinite-component nature of the order parame¬ 
ter. 

There is another evidence of the cr-model defi¬ 
ciency. If the Vollhardt and Wolffe theory is exact, 
then the formalism of dimensional reqularization is 
initially incompatible with the physical essence of 
the problem [2] . The enforced application of this for¬ 
malism should lead to unsolvable problems, and the 
” high-gradient catastrophe” is a possible manifesta¬ 
tion of them. This catastrophe is probably elimi¬ 
nated in other regularization schemes (see discussion 
of the paper [47] in [48]), but a change of the scheme 
surely modifies many-loop contributions. 

We should stress that discrepancy between a- 
models and the self-consistent theory of localization 
arises just on the four-loop level. There is a chance 
for elimination of this discrepancy in the result of 
indicated modifications of er-models. 

8. Conclusion 


We have shown above that multifractal prop¬ 
erties of wave functions can be obtained from 
self-consistent theory of localization by Vollhardt 
and Wolfle, in spite of the opposite statements 
widespread in literature. The diagrammatic inter¬ 
pretation of results allows to derive all scaling rela¬ 
tions used in numerical experiments. Comparison 
with the latter confirms the tendency revealed in 
preceding papers [SiiSi]: the raw numerical 
data are perfectly compatible with the Vollhardt and 
Wolfle theory, while the opposite statements of the 
original papers are related with ambiguity of inter¬ 
pretation and existence of small parameters of the 
Ginzburg number type. 

Analysis of the first e-approximation of the 2 + e 
theory reveals existence of two qualitative proper¬ 
ties: (i) irrelevance of the correlation length £ in the 
metallic phase as a characteristic length scale, and 
(ii) realization of the maximally symmetric form (56) 
for the n-point correlator of wave functions. Due to 
a qualitative character of these properties they have 
a chance to be exact; then the multifractal spectrum 
is strictly parabolic and determined by the one-loop 
Wegner result. A surprising accuracy of this result in 
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application to the d = 3 and d = 4 cases was repeat¬ 
edly reported in literature, while the observed small 
deviations can be explained by the extremely slow 
convergence to the thermodynamic limit discovered 
in Sec.5. The four-loop contribution to the anoma¬ 
lous dimensions is surely deficient and may disap¬ 
pear in result of a modification of the er-model La- 
grangian, which is necessary for taking into account 
the spatial dispersion of D(lu, q) and elimination the 
gradient catastrophe. Simultaneously, it may lead to 
elimination of other discrepancies between er-models 
and the self-consistent theory, which are present on 
the four-loop level. As noted in Sec.3, validity of the 
self-consistent theory is directly connected with the 
property (i). 

As for the relation of multifractality with a spa¬ 
tial dispersion of the diffusion coefficient D(u,q), 
this question is resolved in the compromise man¬ 
ner. A definition of D(ui,q) is ambiguous and al¬ 
lows the ’’gauge transformation”. A spatial disper¬ 
sion is absent in the ’’natural” gauge (88), while the 
Chalker hypothesis mi takes place in the gauge with 
4>reg{q ) = 0 . The raw data of paper [12] on spread¬ 
ing of the wave packet and the physical experiment 
[351 indicate the absence of a spatial dispersion for 
the observable diffusion constant. 
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